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On  Systems  of  Ordinary,  Kon-llnear  Differential  Equations 
Involvinp;  Periodic  Given  Functions  with  a  Small  Period 

Joel  Franklin 

-'-•   Introduction.   In  this  paper  systems  of  ordinary  differential 
equations 


(1.1)       =  f^(t,x^,. ..,x^)      (v=l,...,n) 

dt 


are  considered  in  which  the  given  functions  f   are  periodic  func- 
tions of  time,  t,  with  a  small  period  e  >  0  .   It  is  shown  that, 
if  e  is  sufficiently  small,  then  the  solutions  of  the  system  (1.1) 
are,  in  general,  well  approximated  over  a  finite  time-interval  by 
appropriate  solutions  of  the  system 

dy 

(1.2)       — ^  =  gv(yi''"'yn^      (v=l,...,n), 
dt 

where  the  functions  g   are  the  time-averages 

€ 

By{jl,'",J^)   ^  J    J  ^^it,y-^,.",y^)dt         (v=l,...,n) 

0 

Numerical  estimates  are  given  which,  it  is  hoped,  may  be  useful 
to  the  applied  mathematician.   The  manner  in  which  these  estimates 
may  be  used  is  illustrated  in  the  last  section  by  a  concrete  appli- 
cation to  a  problem  of  a  type  occurring  in  mathematical  physics. 


2,    Assumptions  .   For  any  vector  z  =  (  z-,  ,  Zp,  . .  . ,  z  )  with  real 
components  we  may  define  the  norm  |z|  as  either 

(2.1)  |z|  =  V-  |zj    or  lz|  =  Jy~-  zl]^^^    . 

Let  either  one  of  these  definitions  be  chosen  and  used  throughout. 
Let  T  and  A  be  fixed  positive  numbers.   Now  let  a  vector-function 
f(t,x)  =  (  f-i  ( t,x,  ,  .  .  .  ,x  )  ,  .  .  .  ,f  ( t,x-,  ,  .  .  .  ,x  )  )  with  real-valued 
components  be  defined  in  the  range  0  <  t  <  T,  |x|  <  A.   We  assume 
that  f(t,x)  is  periodic: 

(2.2)  f(t+e,x)  =  f(t,x), 

where  e  is  a  positive  number  less  than  T.   And  we  assume  that 
f(t,x)  is  bounded: 

(2.3)  |f(t,x) 1  <  M  . 

Furthermore,  f(t,x)  should  be  a  continuous  function  of  t  and  x 
except  for  at  most  a  finite  number  of  values  of  t  independent  of  x. 
Finally,  we  assume  that  f(t,x)  satisfies  a  Lipschitz  condition: 

(2.U)  |f(t,x)  -  f(t,y) I  <  L|x-y|, 

where  L  is  a  constant  >  0.   (Bear  in  mind  that  the  estimates  (2.3) 
and  (2.[|.)  are  required  to  hold  only  in  the  given  range  of  definition 
of  f (t,x) . ) 


Kow  for  |>:|  <  A  we  may  define  the  function 

e 

(2.5)  g(x)  =  ^  J  f(t,x)  dt. 

^   0 

It  follows  from  (2.3)  ani  (2.[j.)  that  g(x)  is  bounded  and  that  g(x) 
satisfies  a  Lipschitz  condition: 


(2.6)  lg(x)l<W,    ls(x)-g(y)l  <  L-  |x-y|, 

where   F'    <  l<\   ffnd   L'    <   I    . 

Finally,  we  define  a  scalar  function  p(t)  as  follows: 

t 

(2.7)  p(t)  =  c(t)  +  L   J  e^^^-''"  )a(T  )dr  , 

0 

where 

(2.8)  G(t)  =:  ^  (L+L')M'e[-|]  +  (M+M-)  (|-  [|]); 

here  the  symbol  [-;]  denotes  the  greatest  integer  not  greater  than 
-  .   V/e  note  that  p(0)  =  0;  that  for  t  >  0  the  function  p(t)  is 
continuous  except  for  t  =  e,  2e,  3e,  ...;  that  p(t)  ^  0  for  t  >  0; 
and  that  for  0  <  t  <  T  the  function  p(t)  is  less  than  some  constant 
independent  of  e.   V/e  shall  discuss  this  function  in  greater  detail 
in  Section  I4.. 

Before  proceeding  to  the  theorem,  we  make  the  familiar  remark 
that  the  differential  system  (1.1)  is  equivalent  to  the  integral 


k' 


equatior 


x(t)  =  x^  +   J  f(T  ,  x(T  ))dT  , 


0 

where  x   is  an  arbitrary  constant  initial  vector. 

However,  it  is  to  be  understood  that  the  contir.uous  solution 

X  =  u(t)  need  not  be  dif ferentiable  at  the  values  of  t  for  which 

f  is  discontinuous. 

^*    The  Theoi  em.   First  vje  prove  two  simple  lemmas. 

Lemma  1.   The  function  p(t)  defined  in  (2.7)  satisfies  the 
integral  equation 

t 

(3.1)  p(t)  =  a(t)  +  L  ^  p(T)dT  , 

0 

where  a(t)  is  defined  in  (2.8). 

Proof:      The   ripht-hand    side   of    (3.1)    is    equal   to 

t  t  t 

(3.2)  c(t)    f   L     (    c(T-)dt      +   L^     j  ]    e^-'^'^)a(^)d?dt   . 

0  0   0 

Because  a  is  a  continuous  function  except  for  simple  jump- 
discontinuities,  the  order  of  integration  in  (3*2)  may  be  reversed, 
so  that  the  last  term  in  (3.2)  equals 


t  t 


l2  j  jeL(T-0^(^)^^a^  =  L  J  eL(*-^)a(Od^  -  L 
0  ^  0 


a{K)<iZ      . 


Hence,  the  expression  {3.2)    is  equal  to 

t 
a(t)  +  L  J  e^^^-^\{^)dK   =  p(t)  , 
0 

as  required. 

Lemma  2.   Let  (T  ( t )  be  any  integrable  fxinction  satisfying  the 
Inequality 

t 
(3.3)   cJ-lt)  <  G(t)  +  L  J  c^(T)d       (0  5  t  <  T  )  . 

0 


Then 


(3.1+)       cy(t)  <  p(t)         (0  <  t  <  Tq)  . 

Proof:   Subtracting  (3.3)  from  (3.1),  we  find 

t  V 

(3.5)   p(t)  -  c?(t)  ^  L   \      p(T)  -  cT(T)]  dT   (0  <  t  <  Tq)  . 

0  '• 

Since  the  right-hand  side  is  continuous,  the  function  p(t)  -  (5'(t) 
has  a  finite  lower  bound,  say  -ix  .   Hence,  using  {3.S)    repeatedly, 
we  find 


p(t)  -  (?(t)  ^  -  ^j  L^^Lt^      (k=l,2,3,  ...). 


Letting  k  — >cd,  we  find  p(t)  -  d{t)    ^   0,  as  required. 

Theorem.   Assume  the  conditions  of  Section  2.   Let  x^  be  a 


constant  vector  with  Ix  I  <  A  .   Assume  that  for  0  <  t  <  T  the 

'  o '  -   = 

equation 

t 

(3.6)  y(t)  =  x^  +  j  g{y{r))dT 

b 

has  the  solution  y  =  v(t),  and  that  for  some  a  <  A  we  have 

(3.7)  |v(t)  1  <  a  (0  <  t  <  T)  . 

Let  T^*^  be  the  largest  value  of  t '  <  T  for  which  the  inequality 


max      , ,  >  ^  A-a 

(3.8)         o<t<t'     p^^'  =  -r 


holds.   (Thus,  0  <  T"  <  T.)   Then  the  equation 

t 

(3.9)  x(t)  =  x^  +  5  f(T,x(r))dT 

0 

has  a  solution  x  =  u(t)  for  0  <  t  <  t'"",  and  u(t)  satisfies  the 
inequality 

(3.10)  lu(t)  -  v(t)l  <  ep(t)       (0  ^  t  <  t'")  . 

Remark.   The  sense  of  the  theorem  is  as  follows:   If  e  >  0 
is  very  small,  then  -—  will  be  very  large.   Then  we  should  find 
t""  large,  probably  even  as  great  as  T.   Then  the  inequality  (3.10) 
shows  that  v(t)  closely  approximates  u(t)  over  a  rather  large 
range.   But  the  solution  y  =  v(t)  of  the  initial-value  problem 


It  y(t)  =  g(y(t))  ,     y(0)  =  x^ 


is  generally  comparatively  easy  to  discuss,  since  the  function 
g(y)  does  not  depend  explicitly  on  t,  and  since  g(y),  being  the 
time  average  of  f(t,y),  should  be  in  some  sense  a  smoother  function 
than  f(t,y).   These  observations  will  be  illustrated  by  the  example 
in  Section  S- 

Proof  of  the  theorem.   By  the  classical  existence-and-unique- 
ness  theorem  for  systems  of  ordinary  differential  equations,  we 
know  that  the  equat.'.on  (3  •9)  has  a  unique  solution  x  =  u(t)  in  a 
positive  neighborhood  of  t  =  0,  and  that  this  solution  may  be 
continued  as  long  as  t  and  the  vector  x  =  u(t)  do  not  leave  the 
range  0<t<T,  jxj  <A  .   Let  us  say  that  the  solution  u(t)  may 
be  thus  continued  for  0  <  t  <  T   ,  where  0  <  T^  <  T  .   Then  in  this 
interval  we  have 


(3.11)   ^ 


u(t)  -  v(t)  =  f{f(T,u(T))  -  g(v(T))(  dl 

0 
t  t 

=  I  jfd,  u(T))  -  f(r,v(i))jdt  +   {f(T,v(T))  -  g(v(r))}dT  . 

0  0 

According  to  the  Lipschltz  condition  (2.1+) ,  the  first  of  the  two 
integrals  in  (3.11)  is  not  greater  in  norm  than 

t 
(3.12)         L  J  luCf)  -  v(T)lcir  . 

0 

Letting  N  =  [-]  and  5  =  t  -  eN,  we  may  write  the  second  integral 
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in  (3.11)  in  the  form 

^_^    (k+l)e  t 

(^•^^)  y~  5  {f'(T,v(TM-g(v(T))jdT+  J  {f(T,v(T))-g(v(r))jdT  . 

1{^    ke  t-5 

Since  |f|  <  M  and  |gl  <  M',  the  last  integral  in  (3»13)  has  norm 
not  greater  than  (M+M')5.  Now  the  k--  integral  in  (3.13)  may  be 
written  as 

(k+l)e 

(3.1U)      f   {f(T,v((k4)€))-g(v((k4)e))}  dt 

ice 
(k+l)e 

+    J   {f(T,v(T))-f(T,v((k-ti)e))+g(v((k4l)E))-g(v(T))J  dX  . 

ke 

By  the  periodicity  of  f,  the  first  integral  in  (3  •11+)  equals 


I'  f(T,v((k4)G))dT  -  eg(v((k+^)e))   , 


0 

which  equals  zero,  by  the  definition  (2.5)  of  g.   According  to  the 
Lipschltz  conditions  (2,1|),  (2.6),  the  second  integral  in  (3. Ill-) 
has  norm  not  greater  than 

( k+1 ) e 
(3.15)    (L+L-)     [    |v(T)-v((k+i)e) I  dT  . 

ke 

But  since  v  satisfies  the  equation  (3.6),  and  since  |g|  s.   M',  the 

integrand  in  (3.l5)  is  <  N'  lX-(  k-f^)  €  |  ,  so  that  the  expression  (3.15) 

1         2 

is  <  ■p(L+L')M*e   .   In  summary,  we  have  found 


lu(t)-v(t)l  $L^  |u(i)-v(T)  |clT+  ^(L+L' )M'Ne^  +  (M+M' )5  , 


or 


(3.16)    |u(t)-v(t)|  <  eG(t)  +  L  j  |u(T)-v(r)|dT  , 

0 

where  a(t)  is  the  function  defined  in  (2.8).   We  may,  therefore, 
replace  ^(t)  in  Lemma  2  by  e~  |u(t)-v(t)l,  so  that  {3'k)    becomes 


(3.17)      |u(t)-v(t)l  <  e  p(t)      (0  <  t  <  T^)   . 


Now  it  remains  only  to  show  that  we  may  take  T   equal  to  the  number 
T"  defined  in  the  statement  of  the  theorem.   According  to  the  re- 
marks made  at  the  beginning  of  this  proof,  it  is  sufficient  to  show 
that  every  continuous  vector  x(t)  satisfying  the  inequality 

(3.18)  lx(t)  -  v(t)|  <  e  p(t)       (0  <  t  <  t''') 
must  lie  in  the  range 

(3.19)  |x(t)  I  <  A  (0  <  t  <  t'')  . 


But,  by  the  definition  of  T  "**,  we  have 


P(t)  <  -^-^  ((^  ^  t  <  t'")  . 
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since  Ivj  <  a,  it  follows  that  the  inequality  (3.l8)  implies  the 
required  inequality  (3.19).   We  know,  therefore,  that  the  solution 
X  =  u(t)  of  equation  (3.9)  exists  on  the  whole  interval  0  <  t  <  T"* 
Replacing  T   by  t'"'  in  the  inequality  (3.17),  we  see  that  the  proof 
of  the  theorem  is  complete. 

I4..   Estimation  of  p(  t)  and  t'"'.   In  practice  it  is  often  convenient 
to  have  simple  estimates  for  the  function  p(t)  defined  in  (2.7)  and 
for  the  number  T"  defined  in  the  statement  of  the  theorem.   I'e 
shall  find  an  estimate  for  p(t)  which  preserves  and  makes  obvious 
the  most  important  properties  of  that  function.   However,  our 
estimate  for  t"",  though  simple,  does  not  preserve  the  most  impor- 
tant property  of  t'"',  namely  that  t'"'  >  0  .   (Of  course,  t"  >  0 
because  A-a  >  0,  p(0)  =  0,  and  p(t)  is  continuous  for  0  <  t  <  e.) 
Finally  in  this  section,  we  shall  find  an  exact  representation  for 
p(t)  which  may  sometimes  be  useful. 

To  estimate  p(t),  we  first  note  from  (2.8)  that 

(l|.l)  a(T)  <^(L+L')M'T  +  M+M'  . 

using  this  inequality  in  the  integral  in  (2.7),  we  find 

(i|.2)    p(t)  <  G(t)  +  I^  K'  |e^^-(l+Lt)]   +  (M+M')(e"-1)  , 


or 
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(U.3)    p(t)  <  l(L+L')M'e[7]  +  (M+M')(|  -  [|]) 

+  ^^  M'  ^e^^-(l+Lt)^  +  (M+M')(e^*-1)  . 

To  estimate  T"  ,  we  first  find  a  cruder  estimate  for  p(t)  by 
using  the  inequality  (I4-.I)  in  the  inequality  (I4..2).   We  thus  find 

ik.k)  p(t)  <  (I^M'+H+M')e"  -  I^^  n-  . 

A-a 
Therefore,  certainly  p(t)  <  provided  that  the  right-hand  side 

/*.  -a 
of  (I4..I4.)  is  <  -— -  ,  i.e.  provided  that  t  <  T',  where 

t\     c\  rn,-l<^i^   1.1  o„  l;L(A-a)  ■^   eM'  (L+L'  )  2 

(1^.5)    T'  =  j^   ^log  -  +  log  ^  I.UA   ;  5LM'  f  LMTT  J  • 

It  follows  that- 


(U.6) 

t'"'  >  T'  if  T'  <  T;   T^'  +  T  if  T'  >  T 


Finally  we  compute  p(t)  exactly.   For  convenience,  set 


([^.7)     a^  =  ^i^  ,  a^  ---  ^  (L+L')M'e  -  (H+M'), 


so  that 


([^.8)  c(t)  =  c^t  +  a^l^]      , 
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(i|.9)   p(t)  =  G(t)  +  Le"(a^  | 

0 


e"^rdT   + 


^2  5 

0 


e-LT[IjdT  )  . 


The  first  integral  in  (L|-.9)  is  elementary.   By  the  theorem  of  inte- 
gration by  parts  for  Stieltjes  integrals,  the  second  integral  equals 


-Lel|l 


-  -it|i^-"  ^  i 


-Lke   _        Irti^-Lt   ^   1   ^-Le        1-e 


L'e 


L 


0<ke<t 


1-e 


•Le 


Vie    thus  obtain  the  representation 


-Let|) 


,(t)  =  a(t)  +  ^(eLt-d+Lt)]  -  a^  |  [|l  +  e".  ^'^  ^^  ^  j 

V  \,  j^  ■«  y 


or 


([(..10)   p(t)  =  ^  (e^^-l)  -  a^e^^        ^"^   '  " 


G. 

r 


1-e 


Le 


5.   Example.   For  0  <  t  <  2,  we  consider  a  moving  particle  whose 
displacement,  x, ,  is  governed  by  the  non-linear  differential 
equation 


d2. 


(5.1)   — ^  +  p(t)x^  +  isin^  2007it  +  p(t)l  x^  =  0  , 


where  p(t)  is  the  function  defined  by  the  scheme 


(5.2)   p(t) 


V 


0 

1 

0 

-1 


t  =  0 

0  <  ^  <r^ 

1 


t  = 


1+00 


1  <  t  <  J- 


13 

and  by  the  requirement  of  periodicity: 


(5.3)  P(t)  =  p(t  +  ^) 


The  initial  velocity  of  the  particle  is  known  to  be  zero.   Vie  denote 
the  initial  displacement  of  the  particle  by  the  letter  ^  .   Without 
loss  of  generality,  we  may  assume  that  r  >  0,  since,  if  x,  is  a 
solution  of  the  equation  (5«l)j  so  is  -x,  . 

The  problem  is  to  find  an  upper  bound,  £   ,  for  the  initial 
displacement  r   such  that,  if  ^  <  C   ,  then 

iS^k)  lx^(t)l  <  i^    for  0  <  t  <  2  . 

Of  course,  we  know  from  general  considerations  that  such  a  number 
f   exists;  the  practical  problem  is  to  find  a  value  for  £^  which 
is  not  too  small  in  comparison  with  the  number  -j^,  which  appears 
in  the  inequality  (5.1+)  .   We  shall  find,  by  means  of  the  preceding 
theory,  that  it  is  sufficient  to  take 


(5.5)  Co  =  0-08976   . 


We  shall  also  find  a  closer  description  of  the  displacement  x^  and 
the  velocity  x^  .   In  fact,  if  ?  <  f;^  ,  then  for  0  <  t  |  2 


(5.6)   .-  -  h^t^   -  -^   p(t)  <  x^(t)  ^  .^  -^  ^0 


11^ 

(5.7)   -  I  ^^t  -  ^i^  p(t)  <  ;^(t)  <  ^  p(t)   , 


where 


(5.8)  p(t)  <  0.000258  [200t]  +  0.202  (2C0t  -  [200t]) 

+  0.[^86  (e^-°^^  -  (l+1.06t)j+  0.202(e^-°^^-l). 

To  apply  the  theory  to  the  present  case,  we  first  set 


^  ~  "500  '   ^  "  "To   '    ^  ~  ^ 


Setting   Xp    =   X      ,    we   have 

(5.9)         f(t,x)    =    (X2,    -p(t)x^   -     ^sin^    200Tct    +  p(t)j      x| ) 

Taking  the    time-averafre,    we   find 


1      3 
g(y)   =   (^2    ♦    "  "2  ^1^ 


Next,    we   have 


X     =    {r,    0)    . 
o  \     »        ' 


For  this  problem  we  shall  use  the  second  definition  (2.1)  of  the 

— ^ -x — 

z-j^  +  Z2    .   Accordingly  ,  we 
have 


15 


lg(y)  1^  =  72  ■"  M  =  "'''''  (lO-^sin^G  +  j^   10-  cos  ©)  . 


But 


-^sin^e  +  ^  10-^os^e  =  lo-^-io-^cos^e  +  ^  lo'^os^  5  10'^  . 


10   sin  e  +  r  10   cos  e  =  lu  -iu   cos^w  +  |^ 


Hence , 


M'  =  ^0 


To  find  L' ,  we  write 

|g(x)-g(y) |2  <    (X2-  72)^  ^  li^i  -  yi)'(^i-l^iyil^  yi^ 

^  (x2-y2)'  -^^(-i-yi)'CTfej'  ^  l^-y'"  • 

Therefore , 

L'  =  1  . 

Prom  (5 •9)  we  now  find 

lf(t,x)l2  <  4   ^   x2(l.-2x2)2 

<  maxllO-^sin^e  +  IQ-^cos^G  +  10-2(0.0[^+0.000l;)J 

2 

=  0.0l0)+0l^  =  (0.102)  . 


16 
Hence,  we  may  take 

M  =  0.102   . 

Next, 

|f(t,x)-f(t,y)|2  <  (x2-y2)^  +  (x^-y^)2[l+2(x2+lx^yj+y2)j2 

<  (1.06)2|x-yl2   , 

so  that  we  may  set 

L  =  1.06 

In  summary,  we  have  the  following  data: 

(5.10)  e=0.005,   A=0.1,    T=2, 

L  =  1.06,  M  =  0.102,  L'  =  1,  M'  =  0.1  . 

According  to  the  estimate  ([j..6),  to  ensure  that  T"  =  T  =  2, 
it  is  sufficient  that  the  number^  T'  defined  in  (i+.S)  be  ^  2. 
Replacing  T'  by  2  and  =  by  <  in  (I;. 5),  and  solving  the  resulting 
inequality  for  a,  we  see  that  it  is  sufficient  to  have 

(5.11)  a<A  -.^J-e   (1|LM+5LM'+L'M<  )  -  (L+L')M'?. 
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Computation  with  the  data  (^.10)  gives  (5.11)  the  form 

(5.12)  a  <  0.08980   . 


It  now  remains  only  to  find  the  number  f   in  the  range 


(5.13)  ^   ^   '^o   "   C-08980 


such  that,  if  0  <  ^  <  ^   ,  then  the  solution  of  the  initial-value 


problem 


,2 

^  ^1    1   ^ 
(5.11^)  —2--  +  ^  y^  =  0  ,  y^(0)  =  ■%,   y^(0)  =  0 

dt 


will  satisfy  the  condition 

(5.15)      ly(t)|  <  0.08980        (0  <  t  <  2)  . 

But  the  problem  (5»li^-)  is  very  easy  to  discuss.   Since  y-|(0)  >  0, 
there  must  be  a  number  t   >  0  for  which 


y^(t)  >  0  (0  <  t  <  t^  )   . 


It  now  follows  from  {S-lU-)    that 

(5.16)   y^(t)  <  n,  y^(t)  <  0,  Y^{t)    <   r        (0  <  t  <  t^)   . 

Prom  (5.1i|-)  and  from  the  third  inequality  in  (5.16)  it  follows  that 
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(5.17)  y-^(t)  >  -  1^^,  y^(t)  >  -  -^^^t,  y^(t)  ^4-5  ?^t2(0  <  t  <  t^). 
But, since  0  <  ^  <  f  <  0.1,  we  have 
(5.13)  ''-,   -  T   C^t^  >  0  for  0  5  t  <  I  >  20  . 

Vie   may,  therefore,  take  for  t   any  number  <  20  . 

It  now  follows  from  the  inequalities  (5*16)  and  (5.17)  that 

ly(t)  |2  <  E^   +  (-  ^  ^^t)2  <  ^^   +  r^        (0  <  t  <  2)   . 

Since  0  <  ^  <  0.1,  this  implies 


ly(t)  1  <  E,   /1. 0001      (0  <  t  <  2^  . 

Therefore,  to  ensure  that  a  <  C.o693o,  it  is  sufficient  to  take 
0   <   E   <   r    ,  where  f   is  determined  by  the  equation 

f  i/l.OOOl  =  0.089c30  . 
o  ^ 

Computing,  we  find  the  result  (5.5)*' 

To  obtain  the  inequalities  (5.6)  and  (5.7)  it  is  only  necessary 
to  consider  the  function  p(t)  discussed  in  Section  3.   Using  the  data 
(5.10)  in  the  estimate  (U.3),  we  find  the  inequality  (5.8)  for  p(t). 
But,  according  to  (3.10),  we  have  in  our  case 


x(t)  -  y(t)|  <  ^^   p(t)    (0  <  t  <  2). 
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It    follows    that 
(5.1   9) 

yi(t)  -  2w  P^^^  =  ^1^^)  =  ^1^^^  ^  ^50  P^*^  • 

Applying  the  inequalities  (5.17)  to  the  left-hand  sides  of  (5.19), 
and  applying  the  Inequalities  (5.16)  to  the  right-hand  sides,  we 
obtain  the  required  inequalities  (5.6)  and  (5.7)  . 
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